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Abstract. Streams of gas and dust in the inner parsec of the Galactic center form a distinct feature 
known as the Minispiral, which has been studied in radio waveband as well as in the infrared wavebands. 
A large fraction of the Minispiral gas is ionized by radiation of OB stars present in the Nuclear Star 
Cluster (NSC). Based on the inferred mass in the innermost parsec 10® solar masses), over ^ 10® 
- 10"^ neutron stars should move in the sphere of gravitational influence of the SMBH. We estimate 
that a fraction of them propagate through the denser, ionized medium concentrated mainly along the 
three arms of the Minispiral. Based on the properties of the gaseous medium, we discuss different 
interaction regimes of magnetised neutron stars passing through this region. Moreover, we sketch 
expected observational effects of these regimes. The simulation results may be applied to other galactic 
nuclei hosting NSC, where the expected distribution of the interaction regimes is different across 
different galaxy types. 
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1. Introduction 

The Galactic center hosts the supermassive black hole 
(SMBH) observed as the compact radio source Sgr A*, 
which is surrounded by the Nuclear star cluster (NSC) 
and gaseous-dusty structures, such as HII Minispiral 
arms of Sgr A West, supernova remnant Sgr A East, 
molecular clouds, and the Circumnuclear disk mm- 
It is the closest SMBH and hence its environment 
can be studied with the highest resolution among 
galactic nuclei in the radio-, mm-, submm-, infrared, 
and X-ray wavebands [5]. However, despite high- 
resolution multiwavelength studies several processes 
are still not satisfactorily explained, such as the star- 
formation near the SMBH, the feeding and feedback 
of Sgr A*, and the distribution of the magnetic field 
and its interaction with other stellar and non-stellar 
components. 

The observations of the Galactic center region re¬ 
vealed a large population of young massive stars or¬ 
biting the SMBH as close as ~ 0.1 pc |3]. In fact, the 
NSC seems to be one of the densest concentrations 
of young massive stars in the Galaxy [T]. On the 
other hand, there is an observable flat distribution of 
late-type stars with a radius of as much as 10" [HIS]- 
Thus, a steep relaxed Bahcall-Wolf cusp of stars with 
a slope of 7/4 or 3/2 [317] is probably absent 0 [H] . 

The estimates of the number of stellar remnants 
that use the power-law initial mass function (IMF) 
(standard Salpeter or top-heavy) combined with the 
mass segregation over the age of the bulge (^ lOGyr) 


lead to a considerable population of stellar black holes 
of the order of ~ 10^ [5HTT]. The same order is ex¬ 
pected for neutron stars based on multiwavelength 
statistical studies [H] . Based on the total X-ray lumi¬ 
nosity of the innermost parsec, m set an upper limit 
on the number of compact remnants (< 40 000). 

Such an abundant population of neutron stars ex¬ 
hibiting strong magnetic fields could be utilized to 
further extend our knowledge about the processes in 
the Galactic center. The observations of neutron stars 
(pulsars as well as X-ray sources) near the SMBH 
would contribute to: 

• our understanding of the star formation processes 
near the Galactic center using the number and the 
age distribution of observed sources, 

• mapping the gravitational potential near the SMBH 
using their period derivatives, 

• constraining the electron density profile in the 
Galactic center using their dispersion measures. 

Despite continuing efforts only very few pulsars 
have been detected in the broader Galactic center 
region. It is thought that the lack of detections is 
due to profound interstellar dispersion and scattering. 
However, there are observational hints that such a 
population is present. m report the discovery of two 
highly dispersed pulsars with the angular separation 
< 0.3° from the Galactic center. m confirm the de¬ 
tection of three pulsars with large dispersion measures 
with an offset of ^ 10'-15' from Sgr A*. There is an 
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Figure 1. The distribution of 2302 pulsars in the 
XY Galactic plane. Four logarithmic spiral arms 
are marked by white lines. Sgr A* and the Sun are 
labelled as white points at position (0,0) kpc and 
(0,7.9)kpc, respectively. The data are taken from 
the ATNF Pulsar Catalogue http: //www. atnf . csiro. 
au/research/pulsar/psrcat EZl. 


indication of the relation between the Sgr A East SNR 
and the “Cannonball” source, which was detected in 
both the X-ray and radio wavebands and appears to 
be a pulsar wind nebula [PWN; [TB] . This shows that 
neutron stars can form directly in the Galactic center 
environment. 

In the innermost parsec only one magnetar PSR 
J1745-2900 was affirmed by independent detections 
[IHH2D] at a projected distance of ^ 2.4 ± 0.3" from 
Sgr A*, having a period of 3.7635537(2) [H]. The 
consequence of the lack of sensitivity to detect more 
neutron stars in the Bulge and the Galactic center is 
shown in Figure where most of the detected sources 
are concentrated relatively close to the Sun. 

In this contribution we aim to clarify the distribu¬ 
tion of the interaction modes of neutron stars in the 
Galactic center. Based on the results we discuss the 
possibility of detecting neutron stars in the innermost 
parsec indirectly, specifically by looking for bow-shock 
structures similar to that of the “Gannonball” [T5] . 

The structure of the paper is as follows. In Section]^ 
we explain the set-up of the model and the methods 
that are employed, including the introduction to a 
simple, analytical theory of interaction modes of neu¬ 
tron stars that is, however, sufficient for our purposes. 
Subsequently, in Section]^ we study the distribution 
of interaction modes in the innermost parsec of the 
Galaxy. In Section we discuss the consequences 
of this distribution and the possibilities of detecting 
a fraction of the population indirectly. Finally, we 
summarize our conclusions in Section |5] 

2 . Set-up of the model and 

METHODS 

In our model we concentrate on the innermost parsec 
of the Galactic center, which lies within the sphere 
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Figure 2. The Minispiral model. Top panel: the 
nomenclature of the main components. Bottom panel: 
Keplerian velocity proHle of the gaseons material ac¬ 
cording to |22l I23| . 


of influence of the supermassive black hole (Sgr A*) 
with the radius rsi: 

rsi ~ 1-7 ^ 106Moj UOOkms-iy ’ 

where M, is the central black hole mass and a is the 
stellar velocity dispersion. 

Since we are studying processes that occur at a 
distance of > 1000 gravitational radii from the SMBH 
we approximate the gravitational field by a Newtonian 
point mass of M, = 4.0 x 10® Mq [H |^. We employ 
a Monte Carlo approach for studying the rate of in¬ 
teractions between neutron stars and gaseous-dusty 
structures. Initially, we generate the orbits of neutron 
stars that do not gravitationally interact with each 
other and other stellar and non-stellar components of 
the NSC except for the SMBH. Hence, the orbital ele¬ 
ments of individual stars do not change in the course 
of the simulation. Although this may seem an over¬ 
simplification, it is sufficient for statistical studies of 
the distribution of interaction modes. 

The thermal HII region of Sgr A West known as the 
Minispiral has been studied in the radio-, mm-, and 
infrared wavelengths [^^24j . It consists of four main 
components: the Northern Arm, the Western Arc, the 
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Eastern Arm, and the Bar [see Figure]^ The 

Minispiral consists of a mixture of ionised and neutral 
gas and dust, with temperature ranging from ^ 100 K 
for dust up to lO'^K for the hot ionised gas phase. 
The inferred electron densities are ~ 10^-10^ cm“^ 

MM- 

In our simplified model the Minispiral is represented 
as a system of spherical clumps along the three arms 
20 for each arm) with the length-scale of each clump 
1-2 arcsec. The density and temperature profiles are 
taken from 122H21], but we also analyse the interaction 
modes for values outside this range. The velocity 
profile is assumed to be roughly Keplerian, consisting 
of three streams as inferred from observations by m 
E5] (see Figure [^. We note that deviations from 
Keplerian motion were also detected, probably due to 
magnetohydrodynamic effects [2Z1I23]. 

2.1. Interaction modes of magnetized 

NEUTRON STARS 

Many observed neutron stars exhibit considerable mag¬ 
netic fields, in the case of pulsars most frequently of 
the order of 10^^ Gauss. The warm ionized gas that 
is also detected close to the Galactic centre, specifi¬ 
cally in the Sgr A West region [e.g. HU and references 
therein] has high electric conductivity that is propor¬ 
tional to its temperature, Ac « lO^Te^^cm”^. There¬ 
fore, this plasma must interact effectively with the 
large magnetic field of neutron stars. Gonsequently, 
we need magnetohydrodynamic (MHD) equations de¬ 
scribing the gas dynamics in the potential of a neutron 
star. In fact, if we are close to the neutron star surface, 
we have to solve a system of relativistic magnetohydro¬ 
dynamic equations (RMHD), which is often difficult 
for real systems. 

In our analysis, we focus on the fundamental char¬ 
acteristics of the interaction of a rotating magnetized 
neutron star with the plasma in the surroundings. 
This interaction consists of two parts: gravitational 
interaction characterized by an accretion rate M of 
the captured medium; and electromagnetic interaction, 
which is described by a magnetic dipolar moment /r 
and by a rotational period of a neutron star P. These 
three basic parameters, M, g,, and P, are further com¬ 
plemented by the mass of neutron star Mns and the 
relative velocity with respect to the medium Vao- la 
fact, neutron stars are a part of the broader class 
of gravimagnetic rotators characterized by mass M, 
angular momentum J = /SI, and the magnetic field, 
which is as the first approximation characterized by 
the dipole moment fi. 

In order to simplify the problem, we consider the 
following assumptions: 

(a) the interaction takes place far from the neutron star 

surface, so a relativistic approach is not necessary, 

(b) the electromagnetic part of the interaction is inde¬ 
pendent of the accretion flux parameters. 


(c)the intrinsic magnetic field of neutron stars is a 
dipole field. 



Figure 3. Sketch of a canonical pulsar. In the lower 
right corner, a cross-section of the neutron star is 
depicted with two basic components: the crust and 
the superfluid neutron interior. Magnetic field lines are 
closed up to the light cylinder. Closer to the magnetic 
pole, held lines open and relativistic particles escape. 

Let us make a few notes concerning assumptions 
(a), (b), and (c). Point (a) means that the interac¬ 
tion of interstellar plasma with the magnetic field or 
a relativistic particle wind of the neutron star takes 
place at r GMns/c^, which is often the case. As¬ 
sumption (b) means that an infalling plasma does 
not distort the intrinsic magnetic field of the neutron 
star significantly. Point (c) is approximately valid due 
to the conservation of the magnetic flux. The ratio 
of the quadrupole q and dipole /i components of the 
field scales with the stellar radius as q/g. (x R during 
the core-collapse, so the field is effectively cleansed of 
higher multipole components. The dipole field is the 
most representative component far away from the sur¬ 
face, where the interaction with the ambient medium 
takes place. The closer to the surface the matter gets, 
the more important the quadrupole and the higher 
multipole components become. 

Given assumptions (a), (b), (c), the problem be¬ 
comes much simpler and many effects of the interac¬ 
tion are thus neglected. However, one still gets basic 
information about the character and the scale of the 
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interaction if fundamental parameters of the neutron 
star and those of the ambient medium are given. 

For a non-rotating neutron star, the harmonic sur¬ 
face current j « sin 9 induces a dipole field outside 
the sphere having the following form: 


Bd = 


2/r sin 6 

j'3 


e. 


pcos0 

^3 


eg , 


( 2 ) 


where Gr and eg are unit vectors, fi denotes the dipole 
moment. The angle 9 is measured from the axis 
perpendicular to the dipole field axis (see the sketch 
in Figure [^. The magnetic field magnitude Bq at the 
poles is twice as big as at the equator, which follows 
from the magnitude relation = p/r^(l-|-3sin^ 0)^/^. 
The magnetic dipole moment is given by the field 
magnitude at the poles Bq and the neutron star radius 
Rns, 9- = ^/“^BqR^o, yielding typical values u « 
1/2(2 X 1012 Oe) X 10i8cm3 = lO^OOecm^Q 

In the case of a rotating dipole, the electric com¬ 
ponent of the electromagnetic field also has to be 
considered. If the rotational axis of a neutron star is 
inclined with respect to the dipole axis by angle a, the 
neutron star emits electromagnetic dipole radiation at 
the rotational frequency . In the frame of this model, 
there exists the distance Ri = c/Vt, below which the 
electromagnetic field is static and the electric compo¬ 
nent is given by if ~ v/cB = {Vtr/c)B = (r/R\)B. 
While approaching the radius of the light cylinder, 
i?i, the electric component becomes comparable to 
the magnetic component. On the crossing of the light 
cylinder, the electromagnetic field is no longer static 
and it becomes a propagating electromagnetic wave. 
The magnetic dipole radiation carries away energy 
that may be written in terms of the magnetic dipole 
9, 


difrad 2 2^4 • 2 /n\ 

— = ^^n»ma, (3) 

Gaseous environment and possible modes of 
accretion. A neutron star is assumed to pass 
through an ideally conducting plasma of density poo, 
temperature T^, and the sound speed Coo at infinity. 
We also consider the relative motion of a star with re¬ 
spect to the surrounding medium, v = ri*. The plasma 
starts falling onto the neutron star due to its attrac¬ 
tion. In the case of neutron stars with a negligible 
magnetic field, the stationary flow of matter, or the 
capture rate, is given by the Bondi-Hoyle-Lyttleton 
relation [e.g..E5]: 


, (2GMns)2 „ 


(4) 


where <5 is a dimensionless factor of the order of unity. 
Defining the capture cross-section ctg = SttRq, with 
i?G denoting the capture radius. 


^ 1 Oe is a unit of the H field in cgs units. The equivalent 
in SI units is 1 A/m, 1 A/m = (47r 10“®) Oe. 


i?G 


2GMns 


v'^ + c 


2 ’ 
OO 


(5) 


we may use a convenient form of the capture rate: 


Me = ctgPooT* . (6) 

In realistic astrophysical problems, quantities Too, 
Poo, and Coo are taken at finite distances, R ^ Rq. 

Possible accretion modes may be divided into three 
distinct groups: 

• t;* <C Cooi spherical accretion, without angular mo¬ 
mentum in accreting matter, 

• ti* 3> Coo; cylindrical accretion ~ axially-symmetric 
accretion, without angular momentum in the ac¬ 
creting matter, 

• the disk accretion, infalling matter has a consider¬ 
able angular momentum. 

The following critical radii can be distinguished in 
agreement with the terminology of m- 


Stopping radius. Neutron stars as gravimagnetic 
rotators attract ionized matter due to gravitational 
forces on the one hand, and they prevent accretion due 
to electromagnetic forces on the other. In the rotating 
dipole model, the electromagnetic field is stationary 
inside the light cylinder, Ri = c/fl, and changes into 
a freely propagating electromagnetic wave beyond it. 

The luminosity of the magnetic dipole radiation 
may be estimated by eq. (§, Am — dT^aci/dt. [2)^ 
found out that near the magnetic axis of the neutron 
star, which is inclined by a small angle with respect to 
the rotational axis, the electric component is directed 
along the magnetic field, E « {Q,r/c)BQ, and it accel¬ 
erates charged particles beyond the light cylinder up 
to relativistic energies, which effectively forms a pulsar 
wind. It is assumed that this wind becomes frozen in 
the surrounding medium and passes its impulse to it. 

The pressure of the wind may be estimated as 
Pej = Tm/(47rr2tiej). However, it may happen that 
the accreted plasma penetrates into the region sur¬ 
rounded by the light cylinder. It is prevented from 
accretion onto the surface by the pressure of the static 
magnetic field. Pm = B'^/(8Tr) = p ?To sum 
up, the pressure of electromagnetic forces is of a dif¬ 
ferent character inside and outside the light cylinder, 
and may be approximated by the following relations: 


Pm 



if r < Pi 
if r > Pi. 


(7) 


Using the light cylinder radius. Pi, we may rewrite 
the dipole radiation luminosity of the neutron star 

Tmi 6q. 13 into the form: 




( 8 ) 
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where Kt = 2/3 sin^ a is a dimensionless factor. Using 
eq. ([^, the electromagnetic pressure acting on the 
surrounding gas expressed by 0 may be rewritten 
into the following relations: 


_B. _ 

87rr® 

477 


\i r < R\ 
a r > R\. 


(9) 


Electromagnetic pressure is continuous for Kt = 
1/2 at r = R\. 

Matter being accreted onto the neutron star ex¬ 
erts pressure, which is approximately constant beyond 
the gravitational capture radius Rq and is equal to 
^1‘^Poo L"*- For radii smaller than capture radius Rq, 
matter falls almost freely and exerts dynamic pres¬ 
sure l/2p(r)u(r)^. Under the assumption of spher¬ 
ical accretion, = 47rr^p(r)i;(r), which becomes 
Me = 47ri?QPoo'c* beyond the capture radius, the ac¬ 
cretion pressure Pa can be estimated as follows: 


Pa = 


MeVit 

877 


Me' 




877y>2 \ r J 


if r > Rq 
if r < Rq . 


( 10 ) 


The accretion pressure is a continuous function ac¬ 


cording to eq. (lOl. 


Both the electromagnetic pressure Pm and the accre¬ 
tion pressure Pa as functions of distance are plotted in 
Figure]^ for the following set of parameters: P = Is, 
Bq = 2 X 10^2 G, Rq = 10® cm, nn = lO^cm"®, 
Te = 6 X 10®K, Mns = 1.4Mq, u* = 100kms"®. 

The plasma being accreted by a neutron star is 
halted at the stopping radius where the pressure of the 
electromagnetic forces is in balance with the pressure 
of the accreted matter. The stopping radius may be 
derived from the following condition: 


P — P 

J- m — J- 


( 11 ) 


If the pressure of electromagnetic forces is repre¬ 
sented by the pressure of a static dipole field, then 
eq. (11) yields the Alfven radius. Pa- H the elec¬ 
tromagnetic pressure Pm is related to the outflow of 
relativistic particles (pulsar wind), then the corre¬ 
sponding radius is the so-called Shvartsman radius, 
Psh- The stopping radius may therefore be expressed 
as pH] : 


Pst — 



if Rst < Pi 
if Pst > Pi ■ 


( 12 ) 


For the case when Pgt < Pi, the Alfven radius is 
given by the corresponding relations in eqs. 0 and 
(10), so again we get two cases for distances larger or 
smaller than the capture radius Rq: 


For Pst > Pi, the relativistic pulsar wind interacts 
with the accreted matter. By comparing the depen¬ 
dencies in Figure]^ we see that the accretion pressure 
increases as ~ for Pst < Rq and thus more 

rapidly than the pressure of the pulsar wind, which 
decreases as ^ r“^. Hence, for Pgt < Rq, no stable 
cavern can be maintained by the ejection of matter. 
However, using eq. ([^ for r > Pi and eq. (10) for 
r > Rq, we may find the stopping radius for distances 
Pst > Pg- This radius is also known as the Shvarts¬ 
man radius and may be expressed in the following 
way: 


Rsh = 
Rsh = 


^Tej 

McV*Vej 


1/2 


Rq ) 




MeV^d^ 


(14) 


Let us note that the relations (13) and (14) are valid 


for accretion rates below the Eddington limit for accre¬ 
tion, GMcMNs/Pst < pEdd- The supercritical regime 
occurs in disk accretion modes that do not develop 
when neutron stars are just passing through an ionized 
plasma medium that has low angular momentum. 


Corotation radius. The corotation radius of a ro¬ 
tating neutron star is another important distance scale. 
If accreting plasma penetrates beyond the light cylin¬ 
der, it is stopped at the Alfven radius, Pst ~ Pa, 
where the plasma pressure and the pressure of the 
static magnetic field are in balance. Further evolution 
of stopped plasma is given by the rotational velocity 
of the neutron star. Let us assume that the plasma 
becomes frozen at Pst in the magnetic field and coro¬ 
tates with it at the angular velocity of the neutron 
star, H. The plasma clump will eventually reach the 
surface of the neutron star if the rotational velocity at 
the stopping radius Pst is smaller than the Keplerian 
velocity at the same distance: 


HPst < 




(15) 


However, if the condition (15) is not met, a cen¬ 


trifugal barrier develops that prevents plasma from 
accreting further. The critical corotation radius Pc 
that separates the two modes is given by the equality 


in (15): 


If Pst < Pc, rotation does not considerably effect 
the accretion flux. If Pst > Pc, stationary accretion 
is not allowed. 


Pa 


/ VG^Mg, y/® 

\ MqV^ ) 

( . 

VMc(2GMns)1''V 


2/7 


if Pa > Pg 
if Pa < Rg ■ 


(13) 


Interaction modes. The mode of interaction of 
magnetized rotating neutron stars and also other grav- 
imagnetic rotators is given by the relations between 
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Figure 4. Left panel: the pressure exerted by a stationary magnetic field varies as ~ r“®, whereas outside the 
light cylinder the pressure of relativistic particles that act on infalling plasma falls off with distance as ~ r~^. The 
distance is expressed in light radii (-Ri). Right panel: the pressure caused by accreted plasma is approximately 
constant beyond the gravitational capture radius, whereas for smaller radii gas falls freely and the pressure increases 
as ~ while approaching the star. The units are expressed in gravitational capture radii (Rg)- 


Name 

Notation 

Relation between distances 

Observational effects 

Ejector 

E 

Rst > max{RG,Ri} 

radiopulsars 

Propeller 

P 

Rc < Rst < max{RG,Ri} 

spinning-down more efficient, transient sources 

Accretor 

A 

Rat < Rg and Rst < Rc 

X-ray pulsars, X-ray bursters 


Table 1. Summary of the interaction modes and the types of neutron stars that appear to be relevant in different 
regimes within the environment of the Minispiral. 



Propeller Accretor 



Figure 5. Illustrations of basic interaction modes of magnetized neutron stars with the surrounding environment: 
ejector, propeller, accretor. 


four fundamental distance scales: stopping radius Rgt, 
light cylinder radius Ri, gravitational capture radius 
Rq, and corotation radius Re- Elementary combina¬ 
torics yields 4! = 24 possible relations between these 
radii. However, when taking astrophysics into consid¬ 
eration, we get just a few modes. For example, Ri is 
always greater than Re- We do not take supercritical 
regimes into account in our analysis, since they occur 
in systems with disk accretion, and binary systems are 
also not considered. The classification of three main 
interaction modes - ejector (E), propeller (P), and 
accretor (A) - and thus types of neutron stars accord¬ 
ing to [2^ is summarized in Table The sketches of 
these basic regimes are shown in Figure 


3. Results: Distribution of 

INTERACTION MODES IN THE 

Galactic center 

In our Monte Carlo simulations we produce a suffi¬ 
ciently large synthetic population of neutron stars. 
The number of members of this population is of the 
order of 10^-10®. As a nominal distribution of period 
log P and surface dipole magnetic field log B we take 
the Gaussian fit of the main peak of the observed 
distributions in the catalogue by m 

A(logP) oc exp (^og ^ ^ 

N{logB) oc exp ^ ^ 
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where /ip and /ip are mean values of the Gaussian 
distribution and cTp and ag are the corresponding vari¬ 
ances. The values of the period and surface magnetic 
field are independently attributed to each neutron star 
according to the distribution with the parameters in 
Table E 


Quantity 


(7 

logP 

-0.2188 

0.3488 

logP 

12.0900 

0.4711 


Table 2. Parameters of the synthetic distribution of 
the period (in seconds) and the synthetic distribution 
of the surface magnetic field (in Gauss). 


First, we fix the mean temperature of the ionized 
gas, Te = 6 X 10^ K. The density is continually 
increased by an order of magnitude, starting with 
10^cm“^; we consider densities up to 10®cm“^. 

Neutron stars are distributed uniformly in the space 
of orbital elements, having a uniform distribution of 
the logarithm of the semi-major axis, (0.05,1.00) pc, 
eccentricity, (0.0,1.0), and the cosine of inclination, 
(—1,1). The number of interactions with the Minispi¬ 
ral arms per unit of time is ^ 1% - 10%, depending 
on the length-scale of the individual clumps, 1-2 arc- 
sec, respectively. Thus, for an expected size of the 
population ^ 10^, about 100 to 1000 members may 
interact with the ionized gas of the Minispiral. 

Furthermore, the occurrence of interaction modes 
(E, P, A) is investigated according to the classihcation 
in Table [T] We summarize the results via a series of 
histograms in Figure]^ (Top left panel). It is obvious 
that for typical densities in the Sgr A West region 
inferred from observations [10^-10^ cm“^, refs. [231 HI] 
ejector mode is predominant (> 90%); for larger den¬ 
sities, the fraction of propeller neutron stars increases, 
and at ^ 10^ cm“^ it is approximately the same as 
the fraction of ejectors and starts to prevail for higher 
densities, at which the accretor regime also becomes 
more prominent. This behaviour can be explained by 
the dependence of the stopping radius, eqs. (13) and 
(14), on the accretion rate: for higher densities. Me 
increases and the stopping radius decreases. Thus, we 
naturally obtain a larger fraction of propellers and 
accretors for higher densities. Such densities may hy¬ 
pothetically be found in cold, obscured regions of the 
Minispiral [2S], in the Gircumnuclear disk mm, and 
in dense gaseous-dusty tori and disks of AGN. 

Second, we perform the simulation run with fixed 
density Ug = 5 x lO'^cm”^. The parameters related 
to the magnetic field and the period distributions 
are kept as before. The electron temperature is var¬ 
ied from 1000 K to 60 000 K. We just make a note 
that the motivation for these density and temperature 
variations comes from the fact that the Minispiral 
consists of several streams and clumps of different 
properties that coexist [23]. After several runs, we 
find that the distribution of interaction regimes does 
not vary significantly: ejectors constitute > 92% and 


propellers < 8% of those neutron stars that collide 
with the streams, see Figure]^ (Top right panel). For 
higher temperatures, the accretion rate Mg drops and 
the stopping radius increases, see eqs (13) and (14). 
Hence, the number of propellers decreases whereas 
ejectors are more abundant. However, dependence 
of the distribution on temperature is rather weak. 
No accretors appear, since we keep the density at 
5 X lO'^cm”^ (compare also with Figure]^ - Top left 
panel). 

Moreover, we study the distribution of the inter¬ 
action regimes of neutron stars for stratified Min¬ 
ispiral clumps and different distributions. To this 
end, we implement clumps with a variable density 
profile. They consist of four concentric shells with 
the corresponding outer radii = (i) x 0.08pc, 
N G [1,2,3,4]. The number density increases towards 
the center, starting with 10^cm“^ up to 10^cm“^, 
increasing inwards by one order of magnitude for each 
concentric shell. The temperature is set at the mean 
value of Tg = 6 000K. This scheme ensures that 
most of the gaseous medium has number densities in 
the range (10^, 10®) cm“®, which is in agreement with 
radio and infrared observations. 

We analyse the interaction modes for four different 
distributions whose types and parameters are summa¬ 
rized in Table Distribution (1) is inferred from the 
Gaussian fit of the main peak of both the period and 
the magnetic field distributions of observed neutron 
stars m- Distribution (2) is a broader Gaussian distri¬ 
bution, (3) is a uniform distribution across a large span 
of values: B G (10^ 10^®) G and P G (10-®,102)s. 
Distribution (4) is a combined distribution of (1) and 
(3), where (3) uniformly fills an ellipse in the PB 
plane that satisfies the following condition. 


(logP-/rp)2 {log Bo - ^ 

9 ' 9 _ -L 5 

4 4 

where sp and sp are the corresponding semi-major 
axes, Sp = (logPmax - logPmin)/2 and SB = 
(logs 

max log Benin) 

We perform Monte Garlo simulations with the dis¬ 
tribution of the orbital elements of neutron stars as de¬ 
scribed before. The results of the study of the modes of 
neutron stars that interact with the Minispiral are pre¬ 
sented in Figure(Bottom left panel). In the case of 
distribution (I), the ejector mode dominates, whereas 
the propeller and the accretor modes are represented 
by small amounts, {E; P; A) « (95.4; 4.6; 2 x 10“®)%. 
For the broader Gaussian distribution (2), the ejec¬ 
tor mode decreases, whereas the propeller and the 
accretor modes increase in comparison with the dis¬ 
tribution (1), {E-P]A) « (81.9; 17.8; 0.3)%. This 
trend continues in the case of the uniform distribu¬ 
tion (3), for which the accretor mode rises consid¬ 
erably, {E-,P-,A) « (60.2; 20.3; 19.5)%. Finally, the 
combined distribution (4) with a Gaussian peak re¬ 
sults in the dominant ejector mode again, {E; P; A) « 
(82.3; 13.4; 4.3)%. In all studied cases, the ejector 
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Figure 6. Histograms of the relative distribution of interaction modes of neutron stars in the Galactic center. The 
relative number is in the logarithmic scale. Top left panel: the distribution of interaction modes of the neutron stars 
that interact with the ionized gas in the Minispiral region. Three basic modes are monitored: ejector (E), propeller 
(P), and accretor (A). The distribution alters significantly as the density increases. The temperature is fixed at 
Te = 6 X 10® K. Top right panel: the distribution of interaction modes for the fixed density of fie = 5 x 10"^ cm“®. We 
vary the temperature according to the description in the plot. Only two modes, ejector and propeller, are present. The 
relative number of ejector modes increases (0.9187,0.9189,0.9191,0.9192,0.9192,0.9196) with increasing temperature. 
Bottom left panel: distribution of interaction modes of neutron stars passing through the Minispiral for four different 
distributions: (1) Gaussian, (2) broader Gaussian, (3) uniform distribution, and (4) combined distribution: Gaussian 
-I- uniform. Parameters of individual distributions are summarized in Table Bottom right panel: distribution of 
interaction modes of neutron stars for an increasing distance from the Galactic center. The magnetic field and period 
of neutron stars are distributed according to the combined distribution (4) in Table 


Number 

Distribution type 

Parameters 

1 

Gaussian 

//p = —0.2188, cTp = 0.3488 
fiB = 12.0900, fiB = 0.4711 

2 

Gaussian 

fip = -0.2178, CTp = 0.7019 
= 12.0888, crs = 0.85098 

3 

Uniform 

logPmin = -3, logPmax = 2 
log Hmin = 7, log Hmax = 15 



HP = -0.2189, crp = 0.3478 

4 

Gombined 

= 12.09, o-B = 0.47 
logPmin = -3, logPmax = 2 
log Hmin = 7, log Hmax = 15 


Table 3. Distributions of the period and the maximum surface intensity of the dipole magnetic field of neutron stars 
considered in our work. 


mode dominates (> 60%). However, the distribution 
of the interaction modes is very sensitive to the in¬ 


ternal properties of neutron stars (surface magnetic 
field and period), which evolve with time. The pur- 
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pose of the four discussed distributions of the period 
and magnetic field is to show general trends; more 
detailed studies based on the age distribution of stars 
in the Galactic center are required to constrain the 
distribution further. 

Finally, we study the change in the distribution of 
interaction regimes with increasing distance from the 
Galactic center. We use the combined distribution 
(4) of the period and magnetic field to analyse the 
general behaviour, see Figure(bottom right panel). 
The greater the distance, the smaller the abundance 
of ejectors is, whereas the number of propellers and 
accretors increases. This is due to the fact that the 
relative velocity between stars and gaseous clumps de¬ 
creases for a larger distance, which implies an increase 
in the capture rate Me, see eq. This is translated 
into shrinkage of the stopping sphere, according to 
eqs. (131 and (14). Hence, there is a higher probabil¬ 
ity that the gravitational capture radius and also the 
light radius become larger than the stopping radius, 
which favours the onset of the propeller and accretor 
stage, see Table 


Distribution of interaction modes in the PB 
plane. It is useful to plot the interaction modes of 
neutron stars passing through the ionized Minispiral 
region in the PB plane (period-magnetic field). We 
plot such a map for a uniform distribution (3; see 
Table and distinguish individual regimes by colour: 
ejector - red, propeller - green, and accretor - blue. 
The points in the PB plane where multiple modes 
occur are marked by shades of these colours according 
to the probability of occurrence of particular regimes. 



10 "^ 10 '' 10 “ 10 ^ 10 ^ 
Period [s] 


Figure 7. The uniform distribution of interacting 
neutron stars in the PB plane. Interaction modes 
are color-coded: ejector-red, propeller-green, and 
accretor-blue. Two white lines mark the death lines, 
below which pulsating neutron stars turn off. 

The PB plane is divided by a pair of death lines, 
H « 2.5 X 101° (Q) and H w 3.1 X 10^° (G). 

Above these lines, electron-positron pairs can be pro¬ 
duced in strong magnetic fields, and the pulsars are 
active, whereas below these lines pair production no 


longer proceeds and pulsars are turned off [IS]- The 
ejector mode dominates the region above these lines, 
the accretor stage prevails below them, and the pro¬ 
peller mode occurs in the vicinity of these lines, see 
Figure 

4. Discussion 

Let us briefly discuss plausible observable conse¬ 
quences of the results presented here. 

Ejectors manifest themselves mainly as radiopulsars. 
The predicted number of neutron stars in the central 
parsec of the Galactic centre is ^ lO"!; the expected 
fraction of ejectors is ~ 90% (see Figure]^ assuming 
that the dominant population of neutron stars is the 
same as is observed in the disk population. This 
fraction could be decreased if the population of old, 
isolated accreting neutron stars is more prominent. 
One should be aware of the fact that pulsars stay 
active for a certain amount of time, which can be 
approximated as r « 100 Myr, when they reach the 
death lines and turn off, see Fig. Around the death 
lines they become propellers, and with further spin- 
down and magnetic field decay they start to accrete 
interstellar matter. 

We take the beaming fraction of the ejectors /b = 
0.2[3n],i.e. pulsars that are expected to beam towards 
us, so potentially ~ 1800 pulsars could be detectable 
with sufficient sensitivity. However, there is an appar¬ 
ent lack of pulsar detections in the Galactic centre 
region. The pulsar searches in the Galactic disk are 
generally affected by high background temperature, 
which increases the minimum detectable flux. More¬ 
over, all observations pointed toward the Galactic 
centre region suffer from interstellar scattering result¬ 
ing from turbulent plasma, which causes temporal 
broadening of the pulses to ^ 2000^(7^^ s |3I] (uqhz is 
the observing frequency in GHz). At usual observing 
frequencies ~ 1 GHz it is not possible to detect even 
long-period pulsars, and the surveys have to increase 
the frequency to ~ 10 GHz. However, since pulsars 
have a spectral energy distribution in the power-law 
form Si, cx where a < 0, the flux decreases for 
higher frequencies. To sum up, the lack of observed 
radiopulsars in the Galactic centre may be explained 
by the low sensitivity of current observing facilities. 
A major breakthrough is expected to come with the 
project of Square Kilometre Array-SKA project [32], 
which will provide much higher sensitivity than cur¬ 
rent instruments. Mapping the close vicinity of Sgr A* 
will be the task of the Event Horizon Telescope-EHT 
[55] . which may detect pulsars strongly bound to the 
supermassive black hole, and the data could then be 
used to put the theory of general relativity to very 
precise tests. 

Another way to detect ejectors-pulsars is to look 
for the cavities, pulsar wind nebulae (PWN), and 
bow-shock structures formed when supersonic ejecting 
neutron stars emitting relativistic wind propagate 
through the denser medium of the Minispiral. For 
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Figure 8. The relative velocity of interacting neutron 
stars with respect to the Minispiral streams versus the 
semi-major axes of neutron stars. Red points label 
the interactions with the Northern Arm, green points 
stand for the Eastern Arm, and blue points represent 
passages through the Western Arc. 


the average electron temperature of the plasma in the 
Minispiral region, Te ~ 6 x 10^ K [53], the sound speed 
is Cs « (/cTe/mn)^^^ ~ 7kms“^. According to Fig.[^ 
of the semi-major axis-relative velocity plot, typical 
relative velocities of neutron stars that encounter the 
high-density gas of the Minispiral are greater than the 
sound speed of the interstellar medium by one-two 
orders of magnitude. Hence, the motion is generally 
supersonic, and bow-shock structures associated with 
neutron stars potentially form in the Minispiral arms. 
The size of such a cavity is determined by the energy 
input, E « —E^ot = —/HH, from the central ejector, 
and also by the structure and the properties of the 
surrounding medium. The velocity of the pulsar v^ei 
with respect to the ambient medium also influences the 
size and the orientation of the bow-shock structure. 
Along the relative velocity v^ei the pressure of the 
wind of the ejector and the ram pressure of the ISM 
are in balance at the contact discontinuity. 

The distance of the contact discontinuity from the 
central ejector is given by the stand-off distance [e.g., 

|31j: 


PWN for an ensemble of 10 000 neutron stars, using 
eq. (18) and the assumption that the Gaussian distri¬ 
bution (1) in Table is approximately valid for the 
Galactic center population. Assuming that > 90% of 
interacting pulsars are ejectors, about one thousand 
(in our exemplary runs ~ 1500) bow-shock structures 
could be present along the Minispiral arms. However, 
most of them are rather small (< Imas), and are 
therefore below the detection limit. Propeller bow 
shocks are all below ~ 0.01 mas, whereas the bow- 
shock structures of ejectors have sizes above this value 
and up to ^ 10 mas, see Fig. (Left panel) for compar¬ 
ison; the values along the horizontal axis are decadic 
logarithms of stand-off radii in milliarcseconds. 

The number of larger bow-shock structures (> 
1 mas) is ^ 30 for ^ 1500 neutron stars that inter¬ 
act with the Minispiral out of ^ 10 000 neutron stars 
that occupy the innermost parsec. However, only 
^ 10 of them are greater than 10 mas. Fig. (right 
panel) compares their sizes with the whole Minispiral 
structure, as seen from the Earth for an exemplary 
simulation. 

In addition, neutron stars are expected to occur 
closer to the SMBH, where relativistic effects could 
be detected. They may be members of the S-cluster, 
where they were originally locked in binary systems 
with hot stars. Glose pericenter passages can, however, 
disrupt such binary systems due to triple scattering. 
The components become separated, and for a certain 
set of orbital parameters, one component can escape 
the system on a hyperbolic orbit [see [Ml for details]. 
These hypervelocity neutron stars could potentially 
be seen to interact with the gaseous environment of 
Sgr A West and further away. 

The HII region of Sgr A West, especially the dense 
Minispiral arms and the bar, seem to be a promising 
target to search for the effects of interaction with pass¬ 
ing compact objects, both for the number of encoun¬ 
ters and for the observable effects. Due to its large size, 
the Minispiral appears to be a more favourable target 
for detecting the manifestation of the neutron star 
population than an isolated clump, such as G2/DSO 
currently observed in the Galactic center m- 


rhs 


= ( 


E 


VdTTCPaT^ 


rel 


1/2 


(18) 


where the density of the ambient medium (the Min¬ 
ispiral) is Pa ~ and we adopt the values from 

[531 [53] : the major part of the Minispiral gaseous 
clumps has uh ~ 10'^ - 10®cm“^, with the core of 
each clump of 10^cm“^. Equation (181 describes the 


characteristic size of both the bow-shock associated 
with a fast neutron star and the pulsar wind nebula, 
i.e. the distance of the contact discontinuity from the 
central pulsar. In the case of PWN, the termination 
shock is characterized by concentric tori and jets, as 
is observed in the Crab nebula, for example. 

We compute the typical sizes of bow shocks and/or 


5. Conclusions 

We have studied the distribution of interaction modes 
of the neutron star population in the Galactic cen¬ 
ter. We have explained that the established mode of 
interaction depends on the intrinsic properties of a 
neutron star (the rotational period and the magnetic 
dipole surface intensity) and the external conditions 
(parameters of the gaseous medium: density and tem¬ 
perature) . 

We focused on the distribution of three main inter¬ 
action regimes of neutron stars: ejector (E), propeller 
(P), and accretor (A). We found that the distribution 
of the modes is strongly dependent on the plasma 
density, while it is only weakly dependent on the tem¬ 
perature in the range ^ (10^,10^)K. This further 
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Figure 9. Exemplary simulated distribution of bow-shock sizes associated with neutron stars according to Gaussian 
distribution (1) in Table Left panel: histogram of bow-shock sizes. The green histogram represents propeller 
bow-shocks with rather small characteristic length scales that are beyond the detection limit. Right panel: comparison 
of bow-shock sizes in the simulated 20” x 20” image of the Minispiral. The ejector bow-shocks (red circles) are 
artificially enlarged by a factor of 100 and propeller bow-shocks are enlarged 100000 times enlarged for clarity. The 
elliptical orbits of selected S-stars are also plotted for illustration I3S]. 


implies that the distribution varies across different 
galaxy types. 

We went on to propose an alternative way to look 
for ejecting neutron stars in the Galactic center. One 
can look for the signatures of the interaction of neu¬ 
tron stars with the dense Minispiral arms. Gommon 
structures associated with the supersonic motion of 
pulsars are bow-shocks. Our analysis shows that the 
number of relatively large structures in this region, 
> 10 mas, can be of the order of ^ 10. However, we 
stress that the results depend on the assumed dis¬ 
tribution of the rotational periods and the surface 
magnetic fields. The comparison of different synthetic 
distributions of bow-shock sizes with the observed 
structures will help to constrain the parameters of 
the Galactic center population of young neutron stars. 
These bow-shocks are expected to be sources of polar¬ 
ized non-thermal radiation, possibly detectable in the 
X-ray, infrared, and radio wavelengths. Several cases 
of such comet-shaped sources are indeed observed in 
this region [ISIEH]. Further analysis is required to 
strengthen or reject their association with propagating 
neutron stars. In our upcoming paper we will further 
explore the observability of these features (Zajacek et 
ah, in preparation). 
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